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The main features of the possible quadrupolar glass state in ortho- 
deuterium and para-hydrogen under high pressure are predicted and 
considered in replica-symmetric approximation in analogy with glassy 
behavior of diluted ortho-hydrogen at low pressures. The quadrupolar 
model with J = 2 is suggested. The orientational order and glass regime 
grow continuously on cooling just as it is in the case of ortho-para 
hydrogen mixtures at zero pressure. 

The structure of high pressure phases of solid hydrogen has been widely 
investigated recently (e.g. [1,2, 3]). At low pressures, the centers of para- 
hydrogen and ortho-deuterium molecules occupy the lattice sites of the hep 
structure. The molecules are in spherically symmetric states (phase I or 
SP). Between 28 Gpa for o-D 2 and 110 Gpa for p-H 2 and -150 Gpa the 
molecules become orientationally ordered in broken- symmetry phase II 
(BS). 

High pressure breaks the rotational symmetry of the J=0 solids. The 
ordering occurs even in the systems of para-hydrogen and ortho-deuterium 
molecules because of mixing of higher order J levels with the ground state 
J = as a result of the increase of the intermolecular interactions at high 
pressures [1,2]. The anisotropic interaction between two molecules is 
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dominated at low pressures [5] by electrostatic quadrupole-quadrupole 
(EQQ) interaction, which plays a major role below 100 Gpa [1, 4]. 
However it is necessary to take into account other anisotropic interactions 
to understand the results of precise experiments (e.g. Raman scattering [1]). 

The purpose of this article is to show the possibility of orientational 
glass state in solid p-H 2 and o-D 2 under pressure. Only the main features of 
this phenomenon will be taken into account. The J = —> J = 2 transitions 
take place as a result of the pressure increasing. In this case the J=2 
molecules may be randomly distributed on a close packed lattice. It has 
been indicated in [2, 6] that rotational disorder under pressure may be 
frozen. The molecules with J=0 are spherically symmetric, have no electric 
quadrupole moment and play the role of dilutant. 

So we can consider the quadrupolar glass state in analogy with glassy 
behavior of diluted o-H 2 and p-D 2 at low temperatures and pressures. In 
this case only the ortho-hydrogen and para-deuterium molecular species are 
orientable. They have orbital angular momentum J = 1. The ordered state is 
characteriezed by a long-range orientational order at high ortho-H 2 and 
para-D 2 concentrations [5]. However, for concentrations less than 
approximately 55%, measurements show no evidence of orientational 
phase transition. Instead, NMR experiments have been interpreted in terms 
of freezing of the the orientational degrees of freedom [7, 5]. This state is 
so called quadrupolar glass. 

The EQQ interaction can bring about orientational glass state at high 
pressure in p-H 2 and o-D 2 . The rough estimation of — > 2 transition 
probability a can be done [4, 6] using quantum mechanical perturbation 
theory because the main anisotropic part ( EQQ ) of the intermolecular 
interaction is small [5]. We have a = 4 x 10" 4 for D 2 and a =6 x 10" 5 for H 2 
at zero pressure. The probability to find a J = 2 molecule increases strongly 
with pressure; a = 0.1 at 40 Gpa for D 2 and 150 Gpa for H 2 and a = 0.4 for 
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D 2 at 150 Gpa. It is possible that an intermediate range of a exists where a 
quadrupolar glass occurs. 

Let us consider the system of particles on lattice sites i, j with the 
truncated EQQ Hamiltonian 

» \ -J QO (1) 

where the quadrupole component Q = ^ 

equivalent operator with the same matrix elements Q ~ [3 J] - j(j + 1)] in the 
space J = const. SpQ = 0. 

Quadrupolar glass freezing in mixtures of ortho- and para-hydrogen has 
been considered [8] with 

Q = 3J 2 Z -2,J z = 1,0-1 and J=\. (2) 

This model describes well the zero pressure experiments on ortho- para 
mixtures even in replica- symmetric approach. 

Now the Hamiltonian (1) will be examined on the condition J = 2 and J z 
= 0,1,-1,2, -2. 

2=3^-6], (3) 

J. are random exchange interactions with Gaussian probability 
distribution 




M 2 



can be replaced by 



The scaling J = J /r- , h = ^ °At ensures as usual sensible 



thermodynamic limit. The multiplier 1/3 in (3) is used for simplicity. 



Using replica method the free energy is obtained in the form 
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Extreme conditions for the free energy give the equations for order 
parameters 



where t = J/kT and averaging is performed with the effective Hamiltonian 
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Free energy and order parameters in the replica-symmetric [9] case 
become 
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Here 9j =t 2 ^-^- and 8 2 =zty[q + m 



J 



Order parameters are: m - quadrupolar order parameter (analog of 
magnetic moment in spin glasses), q - glass order parameter and p - 
auxiliary order parameter. 
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exp(4e i )[-exp(-29 2 )+exp(29 2 )]-exp(-9 2 9 1 ) ] 
exp^ )[exp(- 29 2 ) + 2 exp(29 2 )] + 2 exp(-9 2 9! )J ' 
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The temperature dependence of order parameters obtained from (5) - 
(7) is represented in figures 1-3. There is no trivial solution m = 0, q = at 
finite temperature because SpQ 3 * . The orientational order and glass 
regime grow continuously on cooling just as it is in the case of the 
Hamiltonian (1) - (2) describing the ortho-para hydrogen mixtures at zero 
pressure [8]. The quadrupolar long range order is present for T> even if 
Jo = 0. In the pure case (J = 0, J * o) we have from (5) - (6) q 1 2 = m . 

It is easy to see from (2) that Q 2 = 2 - Q . So the equation for p is not 
independent and p = 2-m. There is no similar expression for Q 2 from (3) 
and for order parameter p here. 



Using the equation (4) for the free energy the heat capacity can be 
written in the form 
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Specific heat as a function of (kT/J) calculated from (8) for the four 
cases of (j /j) is shown in figure 4. The dependence of C v on the 
temperature is smooth as in ortho-para mixtures at zero pressure. 

In conclusion, the main features of the random quadrupolar system (1), 
(3) are considered in the replica- symmetric mean-field approximation. This 
system with J = 2 has not been considered earlier. The possible realization 
of the quadrupolar glass state under high pressure in ortho-deuterium and 
para-hydrogen is predicted. 
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Fig. 1 . Order parameters for the case J /J = . 
Fig.2. Order parameters for the case J /J = 1 . 
Fig.3. Order parameters for the case J /J = 2.5. 

Fig.4. Specific heat as a function of {kTIJ) for the cases 

Jjj = 0;J /j=l;J /j = 1.4; JJJ = 2.5 
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Fig. 4 



